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1. Introduction 

PSR B0540-69, a 50 ms pulsar located in the Large Mag- 
ellanic Cloud (LMC), is one of the most distant pulsar s 
known (d = 49.4 ± 3.36 kpc, iTavlor fc Cordesl Il993h . 
ISeward et al. (1984) discovered the pulsar at the posi- 
tion of supernova remnant SNR 0540-693 by using the 
Einstein X-Ray Observatory. Shortly after its discovery 
a bright synchrotr on wind nebula was confirmed in the 
optical waveband ( Chanan et al.lll984l). The optical pul- 
sation s were soon detected bv lMiddleditch fc Pennvpacker 
(1985) with a mean pulsed magnitude of 22.5. In 
the radio band t he p ulsar is q uite a faint source 
( Manchester et all Il993h . though Ijohnston fc Romani 
( 2003f ) recently reported the discovery of the first gi- 
ant pulses from this pulsar. Since 1984, PSR B0540 



69 has been observed by: HST (|Bovd et al 



Serafimovich et al. 1 120041) . V l ^Serafimovich e^al 



Bepp oSAX dMineo et all Il999h. GIN GA dDeeter et al l 
19991) Chandra X-ray Observat ory dGotthelf fc Wand 




2QQ(A iKaaret et all l2Q0lh . ASCA ([ffiravama et al.1 l20Q2h 



and RXTE (|de Plaa et alJl20o3 l 



PSR B0540-69 is often referred to as an extragalactic 
'twin' of the Crab pulsar (Tab. [J). Both pulsars have sim- 
ilar rotational parameters, characteristic age, spin-down 
luminosity and both are embedded in synchrotron pleri- 
onic nebulae. Here the similarity ends. Pulse profiles and 
spectra of PSR B0540-69 and PSR B0531+21 differ signif- 
icantly. The Crab pulse profile shows a sharp double-peak 
structure, whereas the profile of LMC pulsar consists of a 
single broad peak. However, this broad pulse peak might 
be formed as a superposition of two Gaussian compo nents 



separated of about 0.2 in phase (|de Plaa et al.ll20 03). The 
X-ray spec trum of PSR B0540 -60 differs from thet Crab 
spectrum ( de Plaa et al. 2003). Additionally, unlike the 



Table 1. Properties of PSR B0540-69 and PSR B0531+21 
(Crab pulsar). 



Parameter 


PSR B0540-69 


PSR B0531+21 


P(s) 


0.05035 


0.03308 


P (s s- 1 ) 


4.79 x 10" 13 


4.23 x 10" 13 


n 


2.08 


2.51 


d (kpc) 


49.40 


2.0 


r c (kyr) 


1.67 


1.24 


B s (G) 


4.97 x 10 12 


3.78 x 10 12 


L sd (erg s -1 ) 


1.5 x 10 38 


4.6 x 10 38 


Nebula size (pc) 


0.6 x 0.9 


1.5 x 1.5 



Crab pulsar, PSR B0540-69 is not a 7-ray pulsar. 



in X-rays - the Joint European Monitor for X-rays (JEM- 
X, 3-35 keV, Sec. |4j) and in the optical range - the Optical 
Monitoring Camera (OMC, V, 500-600 nm). 

We gathered more than one million seconds of PSR 
B0540-69 observations with INTEGRAL. The data of one 
INTEGRAL orbit, lasting about three sidereal days, are 
summarised as one revolution. The satellite was pointed 
to the LMC during nine revolutions (Tab. [2]); six of them 
took place in January 2003 and three of them a year later, 
i.e. in January 2004. Because of the dithering of the IN- 
TEGRAL satellite, single observations consist of many 
pointings that last about 30 minutes and are separated by 
slews. Each pointing and slew, or only a part of them if 
they are too long, is called a Science Window (ScW). Sin- 
gle observation usually consists of several ScWs. Detailed 
information about the strategy of scientific observations 
and data analysis can be found in the manuals published 
by INTEGRAL Science Data Centre (ISDC) teanJE For 
the purpose of our project we used the INTEGRAL stan- 
dard software - 'Off-line Science Analysis' (OS A, v5.1) to 
generate maps and source intensities. For the timing anal- 
ysis we used our own IDL procedures. 



2. Observations and Data Analysis 

The International Gamma-Ray Astrophysics Laboratory 
(INTEGRAL) is a 15 keV - 10 MeV gamma-ray obser- 
vatory mission, consists of the spectrometer SPI and the 
imager IBIS (Sec. [3]), with concurrent source monitoring 



Detailed description of the data selection 

We used all available INTEGRAL data of PSR B0540- 
69. Until now the LMC was observed during the following 
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Table 2. Time span of the INTEGRAL LMC observations. A total of 813 ScWs within nine revolution were selected. 
First column gives the INTEGRAL orbit (revolution) number. Second column shows the arbitrary numbers of ScWs 
given in brackets, as well as a total number of ScWs per revolution. Next columns show revolution start and end time 
given in MJD TT and UTC, respectively. 



Rev 


Number of ScWs 


TSTART 


TSTOP 


TSTART 


TSTOP 


0027 


(1-92) 


92 


52641.33 


52643.98 


2003Jan02 at 07:57:54 


2003Jan04 at 23:38:55 


0028 


(93-177) 


85 


52644.33 


52646.97 


2003Jan05 at 08:01:21 


2003Jan07 at 23:29:57 


0029 


(178-271) 


94 


52647.33 


52649.97 


2003Jan08 at 07:55:09 


2003JanlO at 23:20:02 


0033 


(272-367) 


96 


52659.29 


52661.94 


2003Jan20 at 06:58:11 


2003Jan22 at 22:34:56 


0034 


(368-459) 


92 


52662.29 


52664.93 


2003Jan23 at 07:09:39 


2003Jan25 at 22:19:59 


0035 


(460-552) 


93 


52665.27 


52667.94 


2003Jan26 at 06:29:29 


2003Jan28 at 22:34:58 


0150 


(553-631) 


79 


53009.26 


53011.93 


2004Jan05 at 06:14:46 


2004Jan07 at 22:29:59 


0151 


(632-726) 


95 


53012.23 


53014.92 


2004Jan08 at 05:37:33 


2004JanlO at 22:14:57 


0152 


(727-813) 


87 


53015.22 


53017.93 


2004Janll at 05:25:23 


2004Janl3 at 22:21:52 



Table 3. Catalogue of selected sources used for the IBIS/ISGRI and JEM-X analysis. Coordinates are given in 
the increasing order of RA. Additionally the ISGRI (LFLAG) and JEM-X (J_FLAG) detection flags are given; - 
undetected, 1 - detected. However, the flag values are obtained for the standard settings of the data analysis performed 
by ISDC team. Especially the minimum detection sigma parameter could play a role here. The source coordinates 
were taken from the INTEGRAL General Reference Catalogue (v. 24) containing all sources detected by INTEGRAL 
an d being brighter tha n 1 mCrab above 1 keV. Only the pulsar position was changed. We used the coordinates given 
bv lKaaret etatl (|200lh on the basis of the Chandra observation. 



NAME 


SOURCEJD 


a[°] 


S[°] 


I.FLAG 


J_FLAG 


LMC X-2 


J052029.2-715736 


80.12167 


-71.96000 





1 


LMC X-4 


J053249.2-662214 


83.205 


-66.37056 


1 


1 


SN 1987A 


J053530.0-691600 


83.875 


-69.26667 








LMC X-3 


J053856.4-640501 


84.735 


-64.08361 





1 


LMC X-l 


J053938.7-694436 


84.91125 


-69.74333 


1 


1 


PSR B0540-69 


J054007.7-692005 


85.04675 


-69.33194 


1 





EXO 0748-676 


J074833.8-674509 


117.1408 


-67.7525 


1 






nine revolutions: 0027, 0028, 0029, 0033, 0034, 0035, 0150, 
0151 and 0152 (Tab. [2]). It resulted in 813 ScWs operated 
in the pointing mode. The exposure time per ScW is typ- 
ically 2300 sec but occasionally ScWs as long as 9000 sec 
or as short as a few hundred seconds are encountered. The 
(arbitrary) ScW numbers corresponding to different revo- 
lutions are given in Tab. O Afterwards, we made the fol- 
lowing selection. Using idxJincH we requested for: ScW be- 
ing of pointing type, the target (pulsar) located inside the 
fully illuminated field of view (that is 9° and 4? 8 for IBIS 
and JEM-X, respectively), and both instrument modes be- 
ing equal to 41, which indicates normal operations mode. 
We obtained two lists of ScWs fulfilling the terms of our 
criteria. The IBIS list consists of 617 ScWs, and the JEM- 
X one of 409 ScWs. None the less, in both cases the final 
number of used ScWs decreased during the analysis. It 
was caused by different reasons that are discussed in the 
according paragraphs. 

Both, IBIS and JEM-X are coded mask instruments, 
therefore it is not possible to deal with one source at a 
time. Each source contributes to the background of the 

2 This procedure allows to select ScWs fulfilling terms of user 
criteria. 



other sources detected in the field of view (FoV) of the 
instrument. Thus, all sources brighter or comparable with 
the target have to be included in the user (input) cata- 
logue. For the purpose of the analysis we created our own 
catalogue. It is based on the general INTEGRAL cata- 
logue (v. 24) and consists of seven sources visible in the 
IBIS FoV. Detailed information of the selected sources is 
shown in Tab. [3j The IBIS FoV is larger than the JEM-X 
FoV, therefore EXO 0748-676 was not visible by the X-ray 
instrument. Consequently only the first six sources listed 
in Tab. [3] were taken into account during the JEM-X anal- 
ysis. 

3. IBIS/ISGRI data analysis and results 

IBIS is a gamma-ray telescope on-board of the INTE- 
GRAL satellite with powerful diagnostic capabilities of 
fine imaging, source identification and spectral sensitiv- 
ity in both continuum and lines. It is able to localise 
weak sources at low energy to better than a few arcmin- 
utes accuracy. Its energy resolution is 7% at 100 keV and 
9% at 1 MeV. IBIS is a 7-ray imager consisting of two 
simultaneously operating detectors, ISGRI and PICsIT, 
covering the energy range from 20 keV to 10 MeV. The 
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Fig. 1. Left: RXTE PCA 2-10 keV pulse profile of PSR B0540-69 obtained form the data collected between December, 
7 2002 and January, 30 2003. Right: RXTE PCA normalised pulse profiles of PSR B0540-69 obtained form the data 
collected between January, 3 and February, 27 2004. The red (lowest), green (middle), and blue (highest) light curves 
corresponds to 2-10, 10-20, and 20-30 keV energy bands, respectively. The pulse duty cycle changes significantly with 
energy. The background level is unknown, therefore the pulsed fraction can not be measured. [Lucien Kuiper, private 
communication] . 



Table 4. PSR B0540-69 ephemeris. 



Parameter & Value Value 



a 2 ooo 05^40 m ll!221 
W) -69°19'54'/98 

Rev. 27, 28, 29, 33, 34 and 35 
52620-52670 
52625.000000000 
19.779298587915 
-1.87410 
2.51 



Val. range (MJD) 
t (TDB MJD) 
iso (Hz) 

is (10" 10 Hz s _1 ) 
Oh (lO -20 Hz s~ 2 ) 



Rev. 150, 151, 152 



Val. range (MJD) 
to (TDB MJD) 
iso (Hz) 

iso (1(T 10 Hz s- 1 ) 
is'o 10" 21 Hz s~ 2 ) 



53007-53062 

53007.000000000 

19.773116586988(11189) 

-1.87256(11) 

3.413(4.92) 



first layer (ISGRI) is made of Cadmium- Telluride (CdTe) 
solid-state detector and the second (PICsIT) of Cesium- 
Iodide (Csl) scintillator crystals. This configuration en- 
sures a good broad line and continuum sensitivity over 
the wide spectral range covered by IBIS. Our analysis is 
based on the ISGRI single events, where the photons are 
stopped in one pixel of the ISGRI layer. The signal am- 



plitude yields the energy of the incident photon. However, 
above 50 keV the energy is a function of not just the pulse 
height but also the pulse rise time (RISE_TIME param- 
eter, afterwords). In that case both measured values are 
used to determine the energy of the photon. 

Lucien Kuiper (private communication) constructed 
pulsar ephemerides from observations with the Rossi X- 
ray Timing Explorer Mission Proportional Counter Ar- 
ray (RXTE PCA). The RXTE data were taken between 
December, 7 2002 and January, 30 2003, as well as be- 
tween January, 3 and February, 27 2004. The ephemerides 
obtained for both epochs are presented in Tab. 2J re- 
spectively, together with the INTEGRAL revolutions to 
which they are applicable. Note that the second frequency 
derivative for the time span 53007-53062 MJD (Rev. 150- 
152) is not significant. The RXTE PCA 2-10 keV pulse 
profile constructed by using the data collected over De- 
cember 2002 and January 2003 and folded accordingly to 
the pulsar ephemeris is shown in Fig. [IJ left panel. The 
pulse profiles obtained from the data collected between 
January and February 2004 in three energy bands: 2-10, 
10-20, and 20-30 keV are shown in Fig.[TJ right panel. It is 
clear that the pulse duty cycle increases, while the pulsed 
fraction likely decreases with energy. The exact number 

Np 



of the pulsed fraction defined as / = Np 1 ^ Dc can not 
be derived because of the unknown background level. The 
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Fig. 2. PSR B0540-69 light curves obtained from JEM-X (left column) and IBIS (middle and right columns) instru- 
ments. The energy bands are indicated in each panel, two cycles are shown for clarity. The PIF values in the plotted 
PIF ranges were multiplied by 10. 



ephemeris for 2003 and 2004 observations are not absolute, 
and the alignment of the corresponding RXTE pulse pro- 
files is not absolute either. Absolute timing would require 
defining the anchor point in the profile and referencing to 
this phase. 

Since October 18, 2004 all public INTEGRAL data 
are available in two formats: revision 1 (Rev.l) and re- 
vision 2 (Rev. 2). In Rev. 2, the correction of all JD time 



stamps for the offsets between the On Board Time (OBT) 
of each instrument is done. In Rev. 2 the data correction 
(COR) step, as well as the instrumental Good Time Inter- 
val (GTI) and deadline handling (DEAD) steps are per- 
formed by the ISDC using OSA (v. 4. 2) at the science win- 
dow level. During the COR step the data are corrected 
for instrumental effects, such as energy and position cor- 
rections, while the GTI step generates, selects and merges 



48 



A. Slowikowska et al.: INTEGRAL observations of PSR B0540-69 



Good Time Intervals to produce a unique GTI that is later 
used for selecting events. The net observing time is also 
computed in this step. Dead and live times are computed 
during the DEAD step of the analysis. The dead time is 
the time during which the instrument was not capable (for 
different reasons) to register the incoming photons within 
a GTI. However, the data correction implemented in OSA 
5.1 is much better. Therefore, it is highly recommended 
by the ISDC team to rerun these three steps (COR, GTI, 
DEAD) over the used data. We performed these correction 
steps over all selected ScWs, which was a very computer 
intensive task. 

It is known from the ISGRI limitations that a prob- 
lem on-board IBIS causes, under some circumstances, 
event times to be shifted by 2 seconds. The OSA soft- 
ware corrects for it, but still the possible jumps might oc- 
cur. Such behaviour can be verified by checking the FITS 
flag keyword TIMECORR. For accurate timing it is rec- 
ommended to use only the data with no possible jump, 
i.e. TIMECORR=0. Therefore, we used only ScWs with 
TIMECORR equal 0. The most up to date version of the 
known limitations, including list of ScWs for which the 
two-second jumps have been detected, can be found on 
the ISDC web pages 0. 

A standard and stand-alone tool ii_pif_built is used 
to calculate the ISGRI Photon Illumination Fraction 
(PIF) for a given source. Once the PIF map is made, 
it is straightforward to filter event lists based on the 
PIF threshold. By applying the PIF filter, one can re- 
duce the background and optimise the signal to noise 
ratio. From experience, we choose PIF> 0.5 as the 
event selection criterion. We run the ii_pif_built proce- 
dure with pulsa r position (a2oqo = 8 5?046754167, £2000 = 
-69?331938889. lKaaret et al.1 (|200lh ) as an input param- 
eter to compute the PIF for each event. By using the 
evts_extract script only single events (evttype=0, 7 < 
RISE.TIME < 90) were selected. Additionally, the noisy 
pixels (SELECT_FLAG=0) were filtered out. The proce- 
dure evts.extract allows to bary centre the events, however 
at this point we did not use it. Firstly, the instrument time 
corrections were applied, i.e. Ill ± 1 0/is and 185 ± 10/is 
for IBIS and JEM-X, respectively ([Walter et al.1 l20Q3h 
and then derived times of arrivals (TOAs) were barycen- 
tred. Afterwards, only events with the energy < 300 keV 
were selected. For these events we stored their TOAs, PIF 
and energy values. TOAs belonging to the first and sec- 
ond set of observations were then folded according to the 
corresponding RXTE ephemeris (Tab. [4]). Hence, we ob- 
tained two sets of the pulsar light curves in five energy 
bands: 17-25, 25-40, 40-100, 100-200, and 200-300 keV. 
The ephemerides are not absolute, therefore to sum up 
the light curves from both observations we used a cross- 
correlation method to determine their relative phase off- 
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Fig. 3. PSR B0540-69 spectrum in 17-300 keV energy 
range. Above 200 keV only upper limit was derived. The 
dashed line indicates the power law fit. 

set. As the template we used the light curve with the high- 
est significance, i.e. the light curve obtained for JEM-X 
data in the 6-10 keV energy band (see Sec.[4j). In the light 
curves only the events with PIF > 0.5 are included, i.e. 
with probability of illuminating the pixel by the source be- 
ing higher than 50%. The results are presented in Fig. [2j 
middle and right columns. 

By standard data processing (ibis_science_ana lysis) we 
produced the sky images for each individual ScW. The 
source is very faint, therefore it was not detected with 
high significance levels in single ScWs. The OSA software 
allows to gather the single ScWs belonging to different 
observations. Thus, we used all selected ScWs to produce 
the mosaic maps. They correspond to the same five energy 
bands used for the light curve extraction (Fig. [2]). In the 
mosaic maps the PSR B0540-69 was detected in the four 
lower energy bands with the following significances: 13.8<j, 
11.2a, 7.6a and 6.2a, respectively (Fig. [4]). There was no 
significant detection above 200 keV. 

For the faint sources the recommended spectral extrac- 
tion method is based on the mosaic sky images. To derive 
the spectrum we used the count rate values obtained from 
these images. We derive the spectrum of the total source 
exclusively from the ISGRI data, treating the calibration 
of JEM-X as not applicable to such a weak source. The 
photon spectrum can be fitted with a power law of in- 
dex - 2.22, which is compatible with the result found by 
( 20061 ). The spectrum and the power law fit 



Gotz et al 



3 http : //isdc . unige . ch/Sof t/download/osa/osa_sw/ 
osa_sw-5 . 1/ osa.issues . txt 



are shown in Fig. [3j 



4. JEM-X data analysis and results 

The Joint European Monitor for X-rays (JEM-X) operates 
simultaneously with the main 7-ray INTEGRAL instru- 
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Fig. 4. The IBIS/ISGRI significance mosaic images of the LMC region {left) and the close-up of the pulsar neighbour- 
hood {right). Maps are for the 17-25, 25-40, 40-100 and 100-200 keV energy ranges from top to bottom, respectively. 
PSR B0540-69 was detected with following significance 13.8a, 11.2a, 7.6a and 6.2a, respectively. 




Fig. 5. The JEM-X significance mosaic images of the LMC region (left) and the close-up of the pulsar neighbourhood 
(right). Maps are for the 3-6, 6-10, 10-15 and 15-35 keV energy range from top to bottom, respectively. PSR B0540-69 
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merits. It is based on the same principles as IBIS and SPI 
instruments: sky imaging accomplished by a coded aper- 
ture mask based on a Hexagonal Uniformly Redundant 
Array (HURA). Its energy range is 3-35 keV, the angu- 
lar resolution at FWHM is 3', the fully illuminated field 
of view is 4.8°, and the time resolution is on the level 
of 122 lis and 1 ms for relative and absolute timing, re- 
spectively. We used the data of the full imaging telemetry 
format only, i.e. the image resolution of the detector was 
256x256 pixels, timing resolution 1/8192 s = 122 /is, and 
spectral resolution 256 Pulse Height Amplitude (PHA) 
channels. The position determination accuracy depends on 
the number of source and background counts and on the 
position in the FoV. The off-axis collimator blocks some 
of the source photons and beyond the fully coded FoV the 
coding is incomplete. There has been a shift in the po- 
sition correction procedure after the beginning of the re- 
processing INTEGRAL revision 1 data to the revision 2 at 
ISDC. Therefore, when using jemx_science_ana lysis we al- 
ways started the analysis from the correction (COR) level, 
even for Rev. 2 data. 

At the beginning 409 ScWs fulfilled our selection crite- 
ria for JEM-X. After the first run of jemx_science_ana lysis 
388 ScWs remained. Some of them showed a 'negative se- 
lection' error, whereas some did not have the fully codded 
imaging events. Similarly to the IBIS/ISGRI analysis we 
used the user catalogue to obtain the images and light 
curves of PSR B0540-69. The list of sources included dur- 
ing the analysis is the same as in the case of ISGRI, except 
EXO 0748-676. This source was not in the FoV of JEM-X 
during the LMC observations. 

The source was not detected above the 3<r level in al- 
most all single ScWs. Distribution of the detection sig- 
nificance of PSR B0540-69.3 for 388 single JEM-X ScWs 
gives: mean 1.35, and standard deviation 0.84. All the indi- 
vidual images from the different science windows gathered 
in the observation group can be combined in the second 
step of image reconstruction, and produce JEM-X mosaic 
images. The combined images have longer exposure time. 
As a consequence, weaker sources which are not visible 
in single ScW can appear in the mosaic images. This is 
the case of the PSR B0540-69, Fig. El The JEM-X light 
curves were obtained in similar way as the IBIS/ISGRI 
light curves. The pulsar profile in the following energy 
bands: 3-6, 6-10, 10-15 keV are shown in Fig. [2] (left col- 
umn). 

5. Summary 

The INTEGRAL analysis presented here is based on the 
INTEGRAL observations of the LMC obtained in Jan- 
uary 2003 and 2004 with a total exposure of ~ 1.5 Ms. 
In the mosaic maps from the total exposure (JEM-X and 
IBIS/ISGRI) a source at the location of PSR B0540-69 is 
clearly visible up to energies of 200 keV. After barycentric 
correction and determination of the pulsar phases, based 



on the ephemeris available from contemporaneous RXTE 
data, the light curves show the characteristic shape of a 
broad pulse up into the 40-100 keV band. At higher ener- 
gies no significant pulsation is detectable. 

For the first time PSR B0540-69 has been detected 
up to ~ 200 keV with pulsations visible up to 100 keV. 
The total source photon spectrum can be fitted with a 
power law with photon index -2.2 and the flux in the range 
17 - 300 keV is about 6 x 10 -11 ergs^cm- 2 . The pulsed 
fraction of the total emission decreases with energy and 
only upper limits could be derived above 1 00 keV assum- 
i ng a l ight curve profile from lower energies. ICheng & Wei 
(|l995h predict in an outer gap emission model for the pa- 
rameters of PSR B0540-69 a significant downturn of the 
synchrotron spectrum around 50 keV, whic h seems to be 
confirmed by our result. Above 100 keV the lCheng fc Wei 
model predicts an inverse Compton spectrum to domi- 
nate, but its intensity would be lower by about a factor of 
10 with respect to the extrapolation from soft X-rays and 
therefore undetectable with the presently available sensi- 
tivity. 
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